function outI = centerForZernike( I, x0,y0, draw )
d = size(I);
outI = zeros(d);
centerX = round(d(1)/2);
siftX = centerX - x0;
centerY = round(d(2)/2);
siftY = centerY - y0;
for j=1:1:d(2)
    yj = j +  siftY;
    if (yj > 0 & yj < d(2))
        
        for i=1:1:d(1)
            xi = i + siftX;
            if(xi > 0 & xi < d(1))
                outI(xi,yj) = I(i,j);
            end;
        end;
        
    end;
end;
if (nargin > 3)
    for i=-1:1:1
        for j =-1:1:1
            outI(centerX + i,centerY + j) = 1.0;
            I(x0 + i,y0 + j) = 1.0;
        end;
    end;
    colormap(gray);
    subplot(2,3,1);
    subplot(2,3,1); imagesc(outI);
    subplot(2,3,2); imagesc(I);
end;

